Spin-orbit-coupling induced domain-wall resistance in diffusive ferromagnets 
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Intrinsically disordered itinerant ferromagnetic alloys are not readily amenable to theoretical 
investigation because of their complexity. The study of disordered magnetic domain walls (DWs) 
is more difficult still because the magnetization direction changes gradually over long length scales. 
In this paper, we investigate diffusive transport through Permalloy DWs, taking into account non- 
coUinearity, alloy disorder, and spin-orbit coupling fully quantum mechanically, from first principles. 
In addition to the effects of magnetization mistracking and anisotropic magnetoresistance, we find a 
new contribution to the resistance of a DW that comes from spin-orbit-coupling mediated spin-fiip 
scattering in a textured diffusive ferromagnet. 
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Introduction. Ferromagnetic alloys such as CoFeB or 
Permalloy (Py), Ni8oFe2o, form the backbone of exist- 
ing magnetoelectronic devices [1] such as spin valves and 
magnetic tunnel junctions [2]. They will play a cen- 
tral role in spin-transfer torque (STT) devices such as 
MRAMS (magnetic random-access memories) [3, 4], spin- 
torque oscillators (STOs) [5-7], and so-called "racetrack 
memories" [8, 9] based upon the "spin-transfer torque" 
(STT) effect whereby a spin-polarized current exerts a 
torque on a magnetization forcing it to precess [5, 10, 11]. 
A realistic description of electrical transport in itinerant 
ferromagnets is made difficult by the degeneracy of the 
partially filled d bands that are responsible for the mag- 
netism and result in complicated Fermi surfaces for or- 
dered materials. The concepts of Bloch states and Fermi 
surfaces that enable the development of transport theo- 
ries for crystals are lost in disordered alloys. Even though 
the spin-orbit coupling (SOC) is in energy terms small, it 
has a large effect on the transport properties of magnetic 
alloys and must be included in any realistic description 
[12]. Finally, the STT effect results when electrical cur- 
rents flow between materials whose magnetizations are 
not collinear; in the case of DWs, the length scales over 
which the magnetization changes is of order 10-100 nm 
[13]. These difficulties all stand in the way of a satisfac- 
tory description of transport properties of Py [14], cur- 
rently the most important candidate for applications. We 
recently extended an efficient scattering formalism of spin 
transport [15] to include SOC that then successfully de- 
scribes the transport properties of alloys such as Py [16]. 
In this Letter, we extend this to treat non-collinearity and 
report on an application to the resistance of Py DWs. 

Most early theoretical studies of DW resistance 
(DWR) focussed on the effect of magnetization mistrack- 
ing, the inability of conduction electrons to adiabatically 
follow an exchange potential rapidly varying in space that 
results in a positive DWR [17]; other mechanisms involv- 
ing impurity scattering in DWs were found to decrease 



the resistance [18, 19]. Much experimental effort has been 
made to identify whether the DWR is positive or negative 
[20-22] and, for the technologically important Py, both 
signs have been reported [23-25]. An additional com- 
plication presented by Py is the anisotropic magnetore- 
sistance (AMR), a dependence of the resistivity on the 
angle between the current and magnetization directions 
that results from SOC. In the early theoretical models, 
SOC was neglected. Recent studies in ballistic metals or 
semiconductors show that SOC gives rise to an intrinsic 
DWR independent of the DW width because the num- 
ber of allowed propagating channels only depends on the 
magnetization direction [26, 27]. In the diffusive regime, 
such a non-local effect is eliminated by disorder scatter- 
ing and how SOC affects the DWR is still unclear. 

The detailed electronic structure of itinerant ferromag- 
nets makes an important contribution to DWR [28] and 
various approaches have been developed to study DWRs 
in particular materials from first-principles [29-32]. Our 
scattering approach takes the full electronic structure 
into account and allows us to treat alloy disorder, SOC 
and non-coUinear magnetism on equal footings and study 
the resistance of Py DWs. A detailed analysis shows that 
three mechanisms contribute to the DWR of diffusive sys- 
tems: magnetization mistracking results in an additional 
resistance inversely proportional to the DW width which 
is only observable in narrow DWs; AMR dominates the 
DWR of wide DWs if there is a component of the mag- 
netization parallel to the current direction that changes 
through the DW; a new adiabatic DWR in diffusive DWs 
that results from SOC-mediated spin-flip scattering and 
is independent of the DW width and profile. It survives 
in the adiabatic limit and is able to distinguish an arbi- 
trarily wide DW from the corresponding collinear ferro- 
magnet. 

Methods. To study electronic transport through a 
DW, we attach semiinfinite (copper) leads to a finite 
thickness of substitutional Ni8oFe2o alloy (Fig. 1) and 
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FIG. 1. (color online). Schematic illustration of the mag- 
netic configurations of (a) Bloch, (b) rotated Neel, and (c) 
Neel DWs. (d) Sketch of the scattering geometry used in the 
calculations in which a finite thickness of Ni8oFe2o substitu- 
tional alloy is sandwiched between semiinfinite copper leads 
and alloy disorder is modelled using 5x5 lateral supercells 
periodically repeated in the x-y plane. Transport is in the z 
direction. 



rotate the local magnetizations inside the scattering re- 
gion to make a 180° DW. We considered three types 
of magnetization profiles, m = [mx{z),my{z),mz{z)], 
corresponding to Bloch [~f{z),~g{z),0], rotated Neel 
[-f{z),0,-giz)], and Neel [g(z), 0, /(z)] DWs. g{z) = 
sech(^^) and f{z) = tanh(^^) for Walker (W) pro- 
files while for linear (L) profiles, f(z) — sin7r( ^^^^ ) and 



cos TT 



Here r^/^ is the DW center and 



9{z) 

^w/L defines the width. 

The electronic structure is first calculated self- 
consistently without SOC for a slab of collinear Py sand- 
wiched between Cu leads using a surface Green's func- 
tion method [33] implemented with a minimal basis of 
tight-binding linearized muffin-tin orbitals (TB-LMTOs) 
[34]. In the scattering region, potentials, charge- and 
spin-densities inside Ni and Fe atomic spheres (AS) 
are obtained using the coherent potential approxima- 
tion. The local density approximation is employed with 
the exchange-correlation functional parameterized by von 
Barth and Hedin [35]. In the transport calculation, the 
self-consistent spin-up and spin-down AS potentials are 
distributed at random in a 5x5 lateral supercell subject 
to the 4:1 Ni-Fe concentration appropriate for Py and 
rotated in spin space as dictated by the DW profile. The 
whole scattering region with a length of 68 nm contains 
some 8300 atoms. The scattering matrix is calculated 
at the Fermi energy using a "wave-function matching" 
scheme [36] also implemented with TB-LMTOs [15] with 
SOC included [16]. The strength of the SOC, Cso is 
determined from the potential gradient within atomic 
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FIG. 2. (color online). DWR of Ni8oFe2o Bloch DWs with 
Walker (a) and Linear (b) profiles as a function of the respec- 
tive width parameters. For each length, we typically consider 
10 different disorder configurations and the error bars are a 
measure of the spread of the results. Insets: DWR replotted 
as a function of 1/\w/l- 



spheres. Within the Landauer-Biittiker formalism, we 
calculate the conductance of the system from the trans- 
mission matrix t as G = {e'^ /h)Ti (tt^) . The DWR is 
defined as Rdw — 1/G — I/Gq, where Go and G are 
the conductances of slabs of collinear Py and of a Py 
DW, respectively. G and Gq are calculated with identi- 
cal random atomic configurations and /c-point sampling 
so as to exclude from i?DW spurious contributions from 
the Sharvin resistances of the leads and from the Cu|Py 
interfaces, and to focus on the resistance arising purely 
from the magnetization rotation. 

Bloch DWs. The DWR is plotted as a function of 
the DW length in Fig. 2 for Bloch DWs whose magne- 
tization rotates in a plane perpendicular to the trans- 
port direction so there is no contribution from the AMR. 
For a narrow DW, incident electrons see a rapidly vary- 
ing (exchange) potential and are reflected by it giving 
rise to a large DWR that increases with decreasing DW 
width. This magnetization mistracking contribution de- 
creases monotonically as the DW width increases and 



vanishes in the wide-DW (adiabatic) hmit in the absence 
of SOC (dashed hnes). As shown in the insets, this con- 
tribution to the DWR is proportional to 1/A for both 
Walker and linear DWs. The additional local resistiv- 
ity due to the DW thus scales with in agreement 
with earlier first-principles calculations [28] and theoret- 
ical models [17, 19]. It was also observed in a recent 
experiment [37] where the DW width in Pt|Co|Pt multi- 
layer was tuned by ion beam irradiation. 

SOC has very little effect on the DWR of narrow Bloch 
DWs which is dominated by mistracking. However, for 
long DWs, the DWR saturates to a finite value of about 
0.1 fi7m^ independent of whether the DW has a Walker 
or linear profile. We will refer to this SOC-related con- 
tribution that survives in the adiabatic limit and has not 
previously been identified as the adiabatic DWR, Ra- 
It distinguishes a DW from the corresponding coUinear 
configuration and contradicts the universal assumption 
that conduction electrons follows a local magnetization 
adiabatically when they flow through a sufficiently wide 
DW. 

The new contribution to the local resistivity is propor- 
tional to the magnetization gradient or 1/A and can be 
understood by generalizing the Levy-Zhang model [17] 
as follows. In the presence of spin texture and SOC, the 
eigenstates are a mixture of spin-up )t) and spin-down 
ID components based on the local quantization axis, i.e. 
[4-+) = a|t) + b\i) and l^--) = -6*|t) + a*\i). If both 
spin texture and SOC are weak, \a\ ^ \b\. A DW leads 
to a contribution to b that is, to leading order in pertur- 
bation theory, proportional to 1/A [17]. In a relaxation 
time approximation, the contribution to the relaxation 
time from disorder scattering can be written in terms of 
a 2x2 local impurity potential v as 
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In the absence of SOC, scattering is spin-conserving so 
v-fi = and spin mixing only results from the non- 
coUinear magnetization, 6 oc 1 /A. As formulated by Levy 
and Zhang, a DW leads to an extra term in the resistiv- 
ity proportional to 1/A^ via spin-conserving scattering 
[17]. SOC makes spin-flip scattering possible, v^^ ^ 0, 
so the leading order correction to l/r depends linearly 
on b. The 1/A^ terms are of higher order and can be 
neglected. Thus the local resistivity in a Py DW has 
the form p(A) — po + Ra/^ + 0(1/A^) and results in a 
constant resistance Ra in the adiabatic limit. 

To conflrm this qualitative picture, we performed cal- 
culations for spin spirals with a fixed pitch, measured in 
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FIG. 3. (color online), (a) Difference between the resis- 
tivities of spin-spiral (p) and collinear (po) Py. The solid 
line is a linear fit to the data with SOC, yielding a slope 
0.102 ± 0.011 fSlm^ in good agreement with the saturated 
DWR in Fig. 2. (b) Calculated saturation value of the DWR 
for Permalloy Walker Bloch DWs as a function of the SOC 
strength. The solid line shows a quadratic fitting. The verti- 
cal dashed line indicates the true SOC strength. 



terms of the rotation At/) of the magnetization between 
adjacent atomic layers, and varying the length of the spin 
spiral. The resistivity is then extracted from a linear fit- 
ting of the total resistance as a function of the length 
[16]. The difference in resistivities of spin-spiral (p) and 
colhnear (po) Py is plotted in Fig. 3(a) as a function of 
A(/), or of the equivalent DW width Xl = 7r/A0. With- 
out SOC, p — Po shows a quadratic dependence on 1/A 
while the relation becomes linear in the presence of SOC. 
A linear fit (red solid line) yields the adiabatic DWR 
Ra = 0.102 ± 0.011 mm^, that agrees wefl with the 
DWR calculations (Fig. 2). 

From Eq. (1) it is not obvious how Ra will depend 
on ^so because of the complicated interplay between the 
non-collinearity and SOC to affect b. To investigate this 
further, we increase ^so artificially. Since the adiabatic 
DWR does not depend on the DW profile, we calculated 
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i?A for Walker Bloch DWs. As shown in Fig. 3(b), the 
adiabatic DWR rises monotonically with ^go exhibiting a 
quadratic dependence up to a value 3.5 times the actual 
value. As the SOC in Py is quite small, this suggests 
that the adiabatic DWRs in materials containing heavy 
elements can be expected to be quite substantial. 

Ncel and rotated Neel DWs. Figure 4(a) shows the 
DWRs of Walker-profile Neel and rotated Neel DWs 
(solid lines) which contain a contribution from AMR. 
For sufficiently small values of Xw, the DWR increases 
with decreasing DW width because of large magnetiza- 
tion mistracking and its behaviour is essentially indepen- 
dent of the DW type. For large values of Xw, the DWR 
decreases for Neel walls and increases for rotated Neel 
walls essentially linearly as a function of Xw- This be- 
haviour can be understood in terms of the AMR. When 
the magnetization is parallel to the current direction, the 
resistivity pn of coUinear Py is 20% larger than p±, the 
value found when the magnetization is perpendicular to 
the current direction. In general, the local resistivity 
depends on the angle 9 between magnetization and cur- 
rent as p{9) = cos^ 9 + p± sin^ 9 [12]. For wide DWs, 
we can estimate the contribution that the AMR makes 
to the DW resistance of Neel (^^amr) rotated Neel 
(^aSr) DWs as {L » A^) 

-Ramr = / dzp[9{z)] - p\\L = -2Xw{p\\ - pi_), 

RImk = / dzp[9{z)] - p^L = 2Xw{p\\ - P±).(2) 

The linear slopes at large Xw in Fig. 4(a) agree with the 
analytical forms in Eq. (2). More quantitatively, we plot 
in Fig. 4(b) the difference between calculated DWRs of 
rotated Neel and Neel DWs as a function of Xw (the solid 
symbols) and they are in perfect agreement with the an- 
alytical form 4Avi/(/9|| — pi_) (the thick line), where 
are obtained from independent calculations for coUinear 
Permalloy [16]. The linear width dependence suggests 
that the DWRs for wide Neel and rotated-Neel DWs are 
dominated by AMR. This picture is consistent with the 
experimentally observed transition from negative to pos- 
itive DWR in Permalloy as the DW width is reduced to 
be atomically narrow [38]. If we subtract the AMR con- 
tribution from the total DWR using the analytical forms 
in Eq. (2), as shown by the dash-dotted lines in Fig. 4(a), 
we find the same width-dependence of DWRs as in Bloch 
DWs-the DWRs decrease with Xw and saturate around 
0.1 iVtw? corresponding to the adiabatic DWR. 

Conclusions. A first-principles study of the resistance 
of Permalloy DWs underlines the importance of SOC. 
The total DWR in diffusive systems originates from mag- 
netization mistracking, an adiabatic DWR, and AMR. In 
narrow DWs, magnetization mistracking leads to a large 
DWR that is independent of the DW type, is inversely 
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FIG. 4. (color online), (a) Calculated domain- wall resistance 
of Ni8oFe2o Neel and rotated Neel DWs with the Walker pro- 
file as a function of the width \w ■ The dash-dotted lines show 
the DWRs after subtracting the AMR contributions described 
by Eq. (2). (b) Difference between the DWRs for Neel and 
rotated Neel walls shown in (a). The thick line shows the ana- 
lytical form 4:\w{p\\ — P±) with p\\{±) taken from independent 
calculations for coUinear Py [16]. 



proportional to the DW width and is dominated by spin- 
conserving scattering. The adiabatic DWR, which arises 
from the spin-Uip scattering as a consequence of SOC, 
does not depends on the DW width. It implies a quali- 
tative difference in the transport properties of a DW in 
the adiabatic limit and a coUinear ferromagnet. In wide 
DWs, the DWR is dominated by AMR which is propor- 
tional to the DW width if there is a change in the angle 
between the current and local magnetization directions. 
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